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- The evolution of the two-point correlation function, ^(r, z), and the pairwise velocity 

Q , dispersion, a{r, z), for both the matter, (^pp, and halo population, (^^h, in three different 

^ I cosmological models: (ilo5-^o)=(liO), (0.2,0) and (0.2,0.8) are described. If the evolution 

I of ^ is parameterized by C{r,z) = (1 + 2;)~('^+^)^(r, 0), where ^(r, 0) = {r/rQ)~"', then 

C^ ; epp ranges from 1.04 ± 0.09 for (1,0) and 0.18 ± 0.12 for (0.2,0), as measured by the 

the evolution of ^pp at 1 Mpc (from z ~ 5 to the present epoch). For halos, e depends 
indeed on their mean overdensity. Halos with a mean overdensity of about 2000 
were used to compute the halo two-point correlation function, ^/j/j, tested with two 



c/3 . different group finding algorithms: the friends of friends and the spherical overdensity 

^ ' algorithm. It is certainly believed that the rate of growth of this ^hh will give a good 

estimate of the evolution of the galaxy two-point correlation function, at least from 
^ ' z ~ 1 to the present epoch. The values we get for Chh range from 1.54 for (1,0) to -0.36 

■ for (0.2,0), as measured by the evolution of from z ~ 1.0 to the present epoch. 

These values could be used to constrain the cosmological scenario. 

The evolution of the pairwise velocity dispersion for the mass and halo distribution 
is measured and compared with the evolution predicted by the Cosmic Virial Theorem 
(CVT). According to the CVT, a{r,z)'^ ~ GQp{z)r'^(,{r, z) or a oc (1 z)'^/'^. The 
values of e measured from our simulated velocities differ from those given by the 
evolution of ^ and the CVT, keeping 7 and Q constant: e = 1.78 it 0.13 for (1,0) or 
e = 1.40 ±0.28 for (0.2,0). 



^This work was done when he was a postdoctoral fellow of the Department of Astronomy, University of Toronto 
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1. Introduction 

The large scale structure of the Universe that we see today is believed to have developed 
from the growth of small perturbations in the matter density driven by gravitational instability. 
The evolution of the clustering of the mass density field depends on the initial conditions via 
the density power spectrum and the mean density of the universe and is therefore a powerful 
constraint on theories of structure formation. The evolution of the galaxy clustering, however, 
need not necessarily follow that of the collisionless component of the mass density field. Galaxies 
have been subject to external phenomena such as tidal interactions, satellite accretion, mergers, 
etc. or internal phenomena such as galactic winds, that it would be unlikely to see the galaxy 
clustering evolution being the same as the clustering evolution of the dark matter. 

Two straightforward statistical tools which describe the clustering properties of galaxies, 
positions and velocities, are the two-point correlation function, ^{x,t), and the pairwise velocity 
dispersion, a{x,t), hereafter comoving coordinates are denoted by x while proper coordinates are 
denoted by r. In a flat universe where initial conditions were generated by a power-law spectrum 
with spectral index n the two-point correlation function should scale as 

ax,t) = as), (1) 

where s = xjf^ and a = 4/[3(3 -|- n)] ([Peebles 19801) . Furthermore, it has been shown that even 
in the case where the hypothesis of scaling is broken, for instance, by a scale-dependent power 



spectrum, relation (1) has proved to be a very good approximation (e.g. Hamilton et al. 1991; 
Padmanabhan et al. 199"^ ). 



In the regime where the density perturbations grow linearly 

e(r,t) = 6(t)2e(r,t,), (2) 

where 6(t) is the growing mode of the density perturbations and £,{r,ti) is the initial correlation 
function {b(ti) = 1). In the particular case of l^o = 1.0, where b is just the expansion factor of the 
universe, a = {l + z)~^, and P{k) oc k"^ : S,{x,z) oc {l + z)^'^x^^"'^^^ or^(r, z) oc (l + 2;)~("'+^V~("+'^). 
On the other hand, in the highly non-linear regime where the hypothesis of stable clustering is 
supposed to work, ^(r, z) = ^(r, 0)(1 -|- A convenient form of parameterizing the evolution of 

e is 

C{r,z) = il + z)--^'+'kir,0), (3) 

as it removes the universal expansion (|Groth fc Peebles 1977| ). If the hypothesis of stable clustering 
is satisfied e = while e = 2 + n m the linear regime (l^o = 1.0 and P{k) oc A;"). 

The present observed galaxy two-point correlation function ^(r, 0) is to a good approximation 



a power-law = ('^/^o) ^ ■ Davis fc Peebles (1983) from the CFA survey find 7 = 1.77 it 0.04 and 



ro = 5.4 lb 0.3h ^Mpc. Loveday et al. (1995) from the Stromlo-APM survey find 7 = 1.71 it 0.05 
and ro = 5.1 it 0.2/i~^Mpc. Values for 7 consistent with those found locally have been measured 



at moderate redshifts (Shepherd et al. 1996; Le Fevre et al. 1996| ) 
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It is seen from the studies by ^hepherd et al. (1996) and Le Fevre et al. (1996) that the 



correlation length, vq, has evolved a great deal. Le Fevre et al. find that vq has decreased by a 
factor of 10, assuming fig = 1, from the present epoch to z ~ 0.6. An €gg ~ 1 it 1 is derived from 
these two studies. 

There exists an extensive literature both observational and theoretical that at least mention 
the name of ^. Here we are interested in those works that deal directly with the time dependency 
of the two-point correlation function for both the density field and the halo population. In 
this sense, the paper by Davis et al. (1985) is pioneer. They studied, among other things, the 
evolution of and '^g^' in different cosmological models, their numerical simulations consisted 
of 32^ particles in a 64^ grid. They found no way of reproducing simultaneously the observed 
amplitude and slope of the correlation function with (^pp in their (r2o,Ao)=(l. 0,0.0) model (a model 
hereafter is represented by a pair of coordinates, where the first coordinate is f^o and the second 
coordinate Aq). A better match was obtained with their low-density models (fiat and open). A 
biased galaxy formation scenario was invoked to save the (1.0,0.0) model and a ^gg was computed. 
Their was always above of their corresponding ^pp. Much of the clustering of these 'halos' 
was due to the pattern imposed by the initial conditions. This work was not intended to study 
the evolution of ^ although it certainly showed it. On the other hand, the evolution they found 
for ^gg is certainly a rough approximation because we know halos do not necessarly arise from 
high peaks and high-peak particles do not necessarly end up in halos. The evolution of has 
also been showed and studied by other authors (e.g. Carlberg 1991 ; Brainerd &: Villumsen 1992| ; 



Brainerd &: Villumsen 1994| ). Carlberg used fof to compute the evolution of ^hh, and found that 
while the correlation length of the density field continues to grow in comoving coordinates, (,hh 
did not change. Brainerd & Villumsen (1994), hereafter BV, with the analysis of simulations of 
128^ particles, in a standard CDM scenario, were deeper in redshift (they started their analysis at 
z = 5 as opposed to z = 2.15 by Carlberg) and found a non-monotonic growth of (,hh- Halos found 
by fof initially produce a biased ^hh (they follow closely the large-scale pattern of filaments and 
sheets just as particles located in high-density peaks did) and then decreases because of mergers. 
The shape of ^hh will be mainly set by four competing phenomena: (1) mergers, (2) formation, 
(3) dynamics, and (4) disruption. By the present epoch, because mergers dominates the scene, 
one would expect to have a (,hh that lies below ^pp, how high will this bias be? is a question whose 
answer will depend on how much are these halos affected by merging and disruption (this will be 
of course an environmental effect). 

In neither of the studies mentioned above was a value for e computed. Recently, two authors, 
at least, have computed the rate of evolution of ^pp using equation (3). Jain (1996) computed e as 
a function of r/ro(a) for a standard CDM scenario. He founds values for e that range from ~ —0.4 
to ~ 2.0. The rate of growth could be much faster that those commonly cited numbers or ~ 1. 
Peacock (1996) found, on the other hand, that the rate of growth in an open universe could solve 
the apparent contradiction that: (1) the slope of seems not to change at all up to z ~ 1 and 
(2) the rate of evolution is relatively rapid with an e value close to 1. 
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In dissipationless N-body simulations a great effort has been made to find collapsed objects 
which subsequently could be asociated with real galaxy halos ( pavis et al. 198^ ; [Carlberg fc 



Couchman 1989| ; Bertschinger & Gelb 1991; Warren et al. 199^ ; |Lacey &: Cole 1994 ; [Summers 



Davis, & Evrard 1995; van Kampen 1995 ; Klypin, Nolthenius, fc Primack 19971 ). In particular, a 
still often used group finding algorithm is friends of friends, fof ( Davis et al. 1985| ). This algorithm 
find groups of particles, halos, that are connected more closely than a specified link length, / (i.e. 
particles that are in an overdensity region in excess of 5min ~ [3"(^/2)^]~^ ~ 2//'^). It is well 
known that it suffers from the defect of joining, once in a while, two halos that are physically 
distinct. To avoid this problem others group finding algorithms were developed. One of these is 
the spherical overdensity algorithm described by Lacey & Cole (1994). Although the identification 
of these collapsed objects with galaxy halos seems at first not a bad approximation, it was soon 



realized that they suffered from the defect of overmerging (e.g. Frenk et al. 1988). Most of these 
halos identified in an early epoch will be destroyed by the time they reach the present epoch and 
this is a function of the mean overdensity of the halo and mass resolution of the simulation (halos 
with the same mass are more tight in simulations with better mass resolution). Some very nice 
algorithms have been invented to solve this problem and here we cite again the paper by Summers 



et al. where a complete discussion about this problem can be found (see also van Kampen 1995 ). 
However, they are very model dependent. In this paper we try to avoid the problem of overmerging 
by considering only halos whose mean overdensity is rather high (~ 2000). We still expect the 
clustering of halos be less than the clustering of galaxies at recent epochs, specially at small scales 
where very massive halos composed of only one particle should contain many galaxies. Therefore, 
the analogy of halos with real galaxy halos should not be taken beyond its scope. The statistics 
of halos or, other galaxy tracers, has been extensively studied at the present epoch and compared 



with observations (e.g. White et al. 1987; Carlberg & Couchman 198£; Gelb & Bertschinger 1994) 



An estimate of the galaxy pairwise velocity dispersion is given by the Cosmic Virial Theorem 
(CVT, Peebles 198q ) 

~ 3J(7)i/(^)2l^o(^)Qro(z)V-T 



a{r, z) 



(4) 



4(7-l)(2-7)(4-7) 

where Q is the three-point parameter (J(1.7) = 4.14). By assuming that Q and 7 do not vary with 
z and that ^ oc (1 + z)~^^^^\ the behavior of a{z) is required to be: cr oc (1 + z)~^/^. An estimate 



of agg{hr = IMpc) from the CFA survey by pavis fc Peebles (1983) gives 300 ± 40km sec 



which in turn produces a value for i7o = Po/Pc {Pc = SHq/SttG = 1.879/i^ x lO^^^g cm~^, h is 
the Hubble constant in units of 100 km sec~^Mpc~^); however, the value of agg can be as high 
as ~ 1000 km sec~^ if Coma is included (|Mo, Jing, &: Borner 1993| ). Because galaxies may not 
dynamically represent the background mass density field, the accuracy of the CVT as an estimator 
of depends on how well galaxies follow the background dark matter. 

This paper is focused to three main goals: (1) do a more systematic analysis of the evolution 
of ^ by computing it for both the density field and a halo population (tested with two group 
finding algorithms), with a high mean overdensity, in three cosmologicals models: (1.0,0), (0.2,0), 
and (0.2,0.8). (2) Include the evolution of the pairwise velocity dispersion. (3) Put results in a 
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convenient manner so as an observer can use them as they are in the paper. The outhne of the 
paper is as fohows. In §2 the characteristics of the simulations are described. In §3 the evolution 
of is discussed and values for the parameter e are given. In §4 the evolution of is discussed. 
In §5 the evolution of the first and second moment of the mass and halo velocity field is presented. 
And finally in §6 a summary is presented. 



2. Numerical Techniques 

The simulations were performed using the adaptive particle-particle, particle-mesh, AP^M, 
N-body code of Couchman (1991). Each simulation consisted of 128'^ particles in a 128'^ grid. 



The initial conditions are generated using the Zeldovich approximation as described by 
Efstathiou et al. (1985). The CDM input density power spectrum, P(A;), is the fit given by 



Efstathiou, Bond, fc White (1992) 



Ah 

Pik) = — ^, (5) 

{1 + [ak + {bk)V^ + icky]''^'' 

where ^ is a normalization constant, a = 6.4/r, b = 3.0/r, c = 1.7/r, and T = 0,Qh. The value of 
h was set to 1 for (1.0,0.0) and (0.2,0.0) while for (0.2,0.8) it was set to 0.7. Therefore T ranges 
from 0.14 for (0.2,0.8) to 1 for (1.0,0.0). The power spectrum is normalized to = 1 in (1.0,0.0) 
and (0.2,0.0), and to COBE in (0.2,0.8). as is the linearly predicted rms mass fluctuation in an 
8h^^ Mpc sphere at the present time. One is left with freedom of choosing any initial value for 
as subject only to the constraint that Zeldovich approximation is applicable. We ran several 
experiments with 64'^ particles. These indicated that by starting a simulation with a rather high 
initial as, as{ti) (a value of 0.3 would be considered high), one might underestimate the present 
pairwise velocity dispersion by up to 30% to 40%, depending on the value of Oq. To make sure 
that our preferred value of 0.1 for as{ti) did not suffer from this effect, we ran two simulations 
with two different values for as{ti): 0.025 and 0.1 for the 0,q = 1.0 model. Negligible differences 
were found in the two-point correlation function and the pairwise velocity dispersion at z = 
between both simulations. All runs are initialized to a{ti = 1) = 1, where ti is the initial time 
in grid units and a is the expansion factor. The force law corresponds to that between two finite 
spherically symetric density clouds with shape given by p{r) = {48/Trr]^){r]/2 - r) for r < r//2, 
where rj is the smoothing length. The value of r/ is constant in proper coordinates, with a value 
at the present epoch, tjq = 50 kpc. The number of timesteps [1,000 for (1.0,0.0), 1,650 for 
(0.2,0.0), and 1,000 for (0.2,0.8)] was chosen sufficiently high to satisfy the stability criteria of the 
numerical integration (see, for example, Efstathiou et al. 1985| ) at all times. In fact, because the 



softening parameter r/ is kept constant in physical units in our simulations, we require only that 
dt ^ min{^J^{ria)^^^ ,3t/2) (expression valid for = 1), where r] and a are in grid units. The 
(0.2,0.0) and (0.2,0.8) both started at as{ti) = 0.1 and had an expansion factor of 27.4 and 14.1, 
respectively. 
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3. Evolution of ipp{r) 

The two-point correlation function was measured by a direct summation of pairs in bins as 
given by the formula 

i{r) = {Np/nN,dV) - 1, (6) 

where Np is the number of pairs between r and r + dr, dV is the volume of this spherical shell, 
is the number of particles taken as centers, and n is the mean number density of particles. The 
correlation function is neither well determined to distances greater than about one tenth of the 
comoving size of the box nor to distances smaller than twice the force resolution. 



3.1. Evolution of and the Parameter e 

The evolution of ^pp in physical coordinates is plotted in Figure || for (1.0,0.0), in Figure |2| 
for (0.2,0.0), and in Figure ^ for (0.2,0.8). The straight line is a power-law of exponent —1.8 and 
epochs are labeled by different symbols: present epoch (+ symbol), z ~ 0.5 (* symbol), z ~ 1.0 
(open circle), z ~ 1.5 (x symbol), z ~ 2.8 (open square), and z ~ 5.0 (open triangle). The symbol 
~ was used because the correct redshift depends on the values of J7o and Aq (see Table H below). 
These epochs have nothing of special. 

If one parameterizes the evolution of ^pp by ^pp{x,z) = [x^^z) / x\^^^'^ (x denotes comoving 
coordinates), where xq{z) is the correlation length as a function of redshift, a value for e can 
be derived by assuming that ro{z) oc (1 -|- z)~^^~^^\ The correlation length and the slope, 7, is 
then computed by fitting a straight line to log^ vs. logx in the range [??, xq]. In Table || are 
shown our results. In the first and second columns are shown the redshift of the model and the 
correlation length in comoving units (as defined by Cpp(xo) = 1). In the third column is shown 
ro in physical units. The correlation length computed by using the fit is shown in the fourth 
column along with its la error, column five. The value of 7 and its uncertainty are located in the 
sixth and seventh columns, respectively. The values of ^pp at 1 Mpc and 0.2 Mpc are shown in 
columns eight and nine, respectively. From this table we see that the correlation length measured 
by assuming that is a power-law is higher than the one defined by ^pp(ro) = 1 (if ^pp were a 
perfect power-law both would coincide), i.e. there appears to be more than one slope (Qq = 1-0), 
one from 1 < ^pp ^ 100 and the second from 100 ^ ^pp < ^r,, where denotes ^pp at the softening. 
(This latter regime likely corresponds to stable clustering.) This effect decreases as we go to earlier 
epochs and this agrees with the idea that the highly non-linear regime bends ^pp. The highly 
non-linear effect decreases as we go to higher redshifts for two reasons: (1) less intrinsic clustering 
and (2) our scheme that fixes the resolution in physical coordinates. 

The e parameter is computed by simply fitting the logro vs. log(l -|- z) (ro as measured by 
the fit), from where one can extract its value by knowing that the slope obtained from the fit is 
just — (e + 3)/7. By taking the mean value of 7 from Table ||, we obtained epp = 1.15 it 0.03 for 
(1.0,0.0), epp = 0.48 ± 0.08 (0.2,0.0), and epp = 0.40 ± 0.22 for (0.2,0.8). The low value of epp for 
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(0.2,0.8) mostly reflects the fact that 7 is rather low at all epochs, in particular the earliest one; 
for instance, epp = 0.96 it 0.20 is obtained when one does not consider the earliest epoch, which 
seems to be a 'transient' epoch. 

A second way to compute e is assuming that in fact we do have a two-point correlation 
function whose shape is independent of epoch and its functionality with 1 + 2; is a power-law. The e 
values we get when the evolution is measured at 1 Mpc, i.e. assuming ^pp(lMpc, z) oc {1 + z)~^^~^^\ 
are: Cpp = 1.04 ±0.09 for (1.0,0.0), Cpp = 0.18 ±0.12 for (0.2,0.0), and Cpp = 0.75 ±0.11 for (0.2,0.8) 
which agree with the above results within la for (1.0,0.0) and within 2a for (0.2,0.0) and (0.2,0.8). 

We have also measured the evolution of (,pp at 200 kpc. If our scenarios are approximately 
self-similar we expect to measure the same rate of growth no matter which scale we use to measure 
it; i.e., we expect to get the same value of e. That seems to apply to (1.0,0.0) model where an e 
value of epp = 1.03 ± 0.14 is computed when evolution is measured at 0.2 Mpc. Notice, however, 
that epp = 0.73 ±0.11 for (0.2,0.0) and epp = 1.21 ±0.17 for (0.2,0.8). This should not be surprising 
since these two scenarios are far from satisfying a simple scaling solution. 

So far, we have not mentioned what to expect for the evolution of ^pp, measured by the 
parameter e, from the theoretical point of view. If the clustering is fixed in comoving coordinates, 
i.e. ^{x,a) = (ro/x)'^, where rg is the correlation length at the present time and r = xa, then 
e = 7 — 3. On the other hand, in a highly non-linear regime, ^ 200, where bound gravitational 
units keep a fixed physical size, the clustering growth is the result of the increasingly diluted 
background and e = 0. In the linear regime, ^ grows as the square of the growing mode of the 
linear density perturbations, b, which in the ^Iq = 1.0 case it is proportional to the expansion 
factor, then ^(x,a) tx a'^x~^ and therefore e = 7 — 1. The e values that we compute from the 
evolution of the correlation length or from the evolution of ^pp at 1 Mpc are slightly higher than 
those ones expected from linear growing (for h ^ a the e value from linear growth is approximately 
the one obtained for = 1-0 multiplied by [5(to)/a(io)]^) where to is the present time): e = 0.94 
for (1.0,0.0), e = 0.11 for (0.2,0.0), and e = 0.7 for (0.2,0.8). 

The evolution of ^pp has been measured from Zg ~ 5 to 2; = 0.0 and the e value would not 
have changed had we measured it with a different Ze value if both self-similarity and equation (3) 
applied. The epp values, along with their error bars, as a function of Ze are shown in Figure ^ 
Our previous results are the particular case of Zg ~ 5. What is clear from this figure is that 
the parameterization given by equation (3) is a very rough approximation of ^(r, z). This is not 
surprising at all since we are measuring the rate of growth in different regimes; for example, 
at 1 Mpc, the epp value increases as Zg decreases, because evolution is passing from the linear 
to non-linear regime, whereas at 0.2 Mpc the epp value decreases because the rate of growth is 
being measured almost only in the highly non-linear regime where we expect e = 0. The higher 
^pp(0.2Mpc, 0) is, the closer one should expect the model be to the stable clustering regime. This 
agrees with what we find for e for the three models. Therefore, e at 0.2 Mpc will be sensitive to 
parameters such as: h, a^ (normalization of P{h)), T. On the other hand, e at 1 Mpc will depend 
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much on the hnear growth of density perturbations which in turn is a function of and Aq- It is 
expected to be much less sensitive to h and fig. 

4. Evolution of ^/i/i 

Halos were initially identified using the friends of friends algorithm. Any other group finding 
algorithm that suffers from the overmerging problem is expected to give similar results. To show 
this we also used for the (0.2,0.8) model the spherical overdensity algorithm, SO, by Lacey &: 
Cole (1994) to identify halos. The mass of each particle is given by nii = ^oPc[ where 
Pc = 2.754 X 10^^ /i^Mo/Mpc^, Lbox is the size of the box in Mpc and N = 128^ is the number 
of particles. Then, for (1.0,0.0) rrii is 1.31 x 10"Mq, while for (0.2,0.0) and (0.2,0.8) they are 
2.62 X 10^°Mo and 1.28 x lO^^M©, respectively. The galaxy-hke mass range used for halos was 
Smj < nifiaioi^Q) ^ 10^^. Because the mass of a particle in a low-density universe scenario is 
lower than in the fig = 1-0 one, the number of particles contained in a halo of a given mass in the 
former scenario is higher by a factor of {h'^Qo)~^- The link length we used was I = 0.1 {I is in 
units of the mean interparticle spacing), which giv6S Si mininiuii niGcin overdensity of (^-mm — 2000, 
for the three models at each epoch. A mean spherical overdensity of 2000 was chosen for the SO 
algorithm. Halos of a rather high mean overdensity were chosen because they presumably suffer 
less merging than their low mean overdensity counterparts. On the other hand, a still higher mean 
overdensity would produce a rather small number of halos, specially at early epochs, and would 
make statistics very uncertain. 

The overdensity, 6i, reached at the time of virialization of a spherical collapse is constant 
for the (1.0,0.0) universe, 178, and it is a function of tcoU for ^ 1 universe, where tcoii is the 
time where the collapse occurs; for example, at ZcoU ~ 5.0, 6i ~ 246 for (0.2,0.0) (it increases with 
time) . An extra set of halos with a variable link length (increasing as we go to earlier epochs) was 
built for the (0.2,0.0) model. The idea was to pick up halos whose difference 5min — at each 
epoch were the same for both (1.0,0.0) and (0.2,0.0). No significant differences were found. 

The evolution of and for the three models is plotted in comoving coordinates in 
Figure |5|, Figure |^, and Figure 0. In Figure we have also plotted the evolution of £,hh as obtained 
by the SO algorithm. The two earliest epochs were drawn out because small numbers of halos 
make ^hh very uncertain. As far as the rate of growth is concerned no significant differences can 
be seen between the S,hh computed by using fof and the one obtained with SO, hereafter results 
are given by using fof as the group finding algorithm. The function ^/j/i(x,z;r2o,Ao) is shaped by 
the interplay of four phenomena: (a) dynamical clustering, the amplitude of ^/^/j increases due to 
gravitational clustering, (b) merging, halos formed at high redshift may merge with other objects, 
increasing their mass out of the selected range, (c) formation, new halos in the selected range are 
formed by accretion or merging, and (d) disruption, some halos are destroyed by tidal forces or 
two-body encounters. Several characteristics can be detected in Figure |5|, Figure ^ and Figure |^ 
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(1) ^hh initially traces the filamentary structure of the universe; i.e., halos are born "naturally" 
highly clustered, producing a biased S,hh- (2) The clustering of halos does not grow continuously, 
contrary to the clustering of the mass density field. (3) If a lower mean overdensity were used we 
would have a greater overall suppression of (,hh- This would reflect the fact that low-overdensity 
regions are more subject to merger and accretion, which decrease ^hh by moving halos out of 
the selected range ( parlberg 1991| ). (4) The evolution of is essentially fixed in comoving 
coordinates. This is specially true for the models of low-density. (5) Dynamical clustering does 
seem to drive the clustering of halos at later epochs at scales x ^ 0.7Mpc, specially in the (1.0,0.0) 
model. 

4.1. Evolution of S^hh and the Parameter e 

A Table § similar to Table |l| was built for ^hh- In this case the slope of ^hh was determined 
using the separation range 1.0 < a;(Mpc) < 10. Unlike Table [l|, in these tables we have preferred 
to show a mean ^/j^. The la error bars are computed as follows: a set of 10 realizations is built 
at each epoch and for each model, each one being a random subset of its corresponding halo 
population. The number of halos in each subset is one third of its corresponding halo set. A 
mean ^/^j^ with its typical deviation is then computed. As was shown in Figure I and Fi gure p the 
comoving correlation length at first decreases, reaches a minimum and then starts increasing. This 
behavior does not seem to depend on Qq (or on the mean overdensity of the halos). 

In Figure ^ we have plotted the evolution of ^^/^ (open squares) at 1 Mpc and 0.2 Mpc. The 
reason why we are choosing 200 kpc is that it is the smallest scale where we could still measure ^hh 
with a certain degree of confidence, i.e. not being subject to statistical noise produced by having 
too few particles. The evolution of is marked by asterisks. 

In Figure ^ we have plotted ehh as a function of z^- This Figure is the corresponding to 
Figure ^ but for halos. It is interesting to see that the evolution at 1 Mpc from z = to Zg ~ 1.3 
coincides with that of the mass density field for (1.0,0.0) model, in agreement with the idea that 
in this redshift range halo clustering is essentially driven by the dynamical clustering of the mass 
density field. This is not sustained at small scales where mergers make the rate of growth of (,hh 
lie below the rate of growth of ^pp. This was not discussed by BV because they did not measure 
(,hh at these small scales. In the low-density models merging still appear to play a big role in 
shaping the evolution of (,hh- At scales of 200 kpc the growth of (,hh is almost fixed in comoving 
coordinates for both scenarios, i.e. e ~ 7 — 3: for (1.0,0.0) 7 ~ 2.2 and ehh ~ —0.8, for (0.2,0.0) 
7 ~ 1.6 and ehh ~ —1-4, and for (0.2,0.8) 7 ~ 1.5 and ehh ^ —1.5. It is clear from what has been 
discussed so far that evolution in the two-point correlation can be a clear discrimation between 
models of low and high Qq. However, little can be said between models with and without Aq 
unless the evolution of ^gg is closer to that of the density field. 
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5. Evolution of oi^pp and (Ji^hh 
The 1-D pairwise velocity dispersion is defined to be 

^?=<h-^ii)'>, (7) 

where v\\ stands for the component along the line connecting the pair of the relative vector 
velocity. Quantities under a bar or in brackets denotes mean. We have measured (Ti^pp(lMpc, z) 
and computed a value for e assuming that CFi^pp oc (1 + The values for e we get are: 

e„ = 1.78 ± 0.13 for (1.0,0.0), = 1.40 ± 0.28 for (0.2,0.0), and = 2.72 ± 0.28 for (0.2,0.8) 
(evolution is measured from z ~ 5 to z = 0.0). They are higher than the e values we get from the 
evolution of the correlation function. Our simulations find more dynamical evolution than that 
predicted by the the CVT under the assumption that the parameters Q and 7 do not change with 
epoch. 

It is interesting to see that our values for the present 1-D velocity dispersion (^ 1000 km sec~^ 
at 1 Mpc) for Q.q = 0.2 are higher than some values previously calculated in the literature (e.g. 
Davis et al. 198^ ; Kauffman &: White 199^ ). This, we believe, is due to three things: (1) the 



higher resolution of our simulations, (2) the (Ts{ti) effect mentioned earlier (see section §2), and (3) 
our rather high normalization. High values for o"i^pp at small scales have been previously measured 
(e.g. Martel 1991). The value for the velocity dispersion obtained by our open model is more than 
three times the observed galaxy velocity dispersion. 

Table ^ and Table ^ show the evolution of the first and second moment of the mass and halo 
velocity field for our models. We have preferred to show mean quantities for the halo population, 
measured as we did for the two-point correlation function. As expected, less evolution is measured 
in cTi fih as compared with cii^ppbecause halos are subject to merging and accretion. On the other 
hand, the picture we get from the biased galaxy formation scenario, where there is a time when 
galaxy formation ceases (or slows down) followed by dynamical clustering evolution, appears to be 
just part of the story. It is here precisely where one expects halo dynamics to trace better galaxy 
dynamics, and this view is supported by the evolution of S^hh- The pairwise velocity dispersion, 
however, decreases from z = 0.5 to z = 0.0 in all models. This may be due to the high values of 
infall velocities encountered at z = 0.5. 



6. Summary 

We have measured the evolution of the two-point correlation function and the pairwise 
velocity dispersion of the mass density field and halo population. The evolution is parameterized 
mostly by the e parameter. Our e values depend on the scale and the time period where evolution 
is measured, and for halos, they also depend on their specified mean overdensity. Results were 
quoted just for a mean overdensity of about 2000 because we believe halos of this overdensity suffer 
less from the overmerging problem. The e values for S,pp range from 0.4 (for (r2o,Ao) = (0.2,0.0)), 
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when evolution is measured at 0.2 Mpc, to 1.5, when evolution is measured at 1 Mpc (for 
(Qo,Ao)=(l. 0,0.0), both covering a period of time from z ~ 1 to present epoch. The range of €hh 
values covered by halos is: -1.0 ^ ehh ^ 1-5 for (1.0,0.0), -1.4 ^ ehh ^ -0.3 for (0.2,0.0), and 
-2.1 ^ ehh ^ -0.6 for (0.2,0.8). The de gree to which these results constrain the mean density of 
the Universe depends on how well the evolution of the galaxy clustering is traced by the evolution 
of the mass density field or halo population. More and better observations of ^gg at diferent 
rcdshifts along with better numerical determinations of its evolution are needed to constrain more 
the cosmological parameter space. 

The correlation length, ro{z), was computed in two ways: (1) using the "standard" definition 
^{ro,z) = 1 and (2) fitting ^ to a power-law function at each epoch; i.e., assuming that 
^(r, z) = (r / rQ{z))^"' , the shape of is a power-law and does not change with time. The range of 
data used for the fit was from the softening length to the correlation length (as found by the first 
method), to correspond to observations. We like to point out that the analytic evolution does not 
really do a good job of predicting the details of the nonlinear correlation function, errors as large 
as 50% are obtained, so these simple fits are considerable valuable. 

The evolution of (7(lMpc) was assumed to be a power-law with an exponent — e^/2 (justified 
by the Cosmic Virial Theorem, CVT). All values found for e„ (halo and mass density field) are 
systematically higher than those predicted by the CVT; i.e. we see more evolution in the velocities 
than that predicted by the evolution of ^ and the CVT for constant Q and 7. 

PC acknowledges a fellowship from Direccion General del Personal Academico, DGPA, 
UNAM, support from NSFRC at the University of Toronto, and a grant from DGPA through the 
project IN-109896. RGC and HMPC were supported by NSERC. We thank Piotr Zembrowski for 
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Table 1. Evolutionary Parameters of 



z 




ro ^ 






7 




ipp (1 Mpc) 


ipp (0.2 Mpc) 












,Ao)=(1.0,0.0) 






0.00 


4.63 


4.63 


4.86 


0.14 


1.86 


0.05 


35.958 


645.54 


0.52 


3.17 


2.09 


2.20 


0.06 


1.94 


0.06 


5.490 


163.33 


0.91 


2.62 


1.37 


1.42 


0.03 


1.98 


0.06 


1.933 


71.51 


1.48 


2.12 


0.86 


0.85 


0.02 


2.02 


0.05 


0.745 


23.84 


2.85 


1.37 


0.36 


0.35 


0.01 


2.12 


0.03 


0.150 


2.85 


4.87 


0.85 


0.15 


0.14 


0.00 


1.74 


0.07 


0.025 


0.64 












,Ao)=(0.2,0.0) 






0.00 


6.05 


6.07 


6.06 


0.12 


1.91 


0.04 


45.48 


1000.1 


0.50 


4.99 


3.32 


3.37 


0.07 


1.91 


0.04 


9.72 


263.8 


n QQ 
u.yy 


/I nfi 




Z. ID 


u.uo 


1.91 


0.05 




lino 
iiu.y 


1.40 


3.68 


1.54 


1.59 


0.03 


1.88 


0.03 


2.00 


52.3 


2.74 


2.47 


0.66 


0.70 


0.02 


1.88 


0.03 


0.56 


8.2 


5.00 


1.49 


0.25 


0.27 


0.00 


1.67 


0.05 


0.14 


1.4 












,Ao)=(0.2,0.8) 






0.00 


6.42 


6.42 


9.31 


0.53 


1.55 


0.06 


44.21 


453.6 


0.49 


5.33 


3.58 


4.03 


0.21 


1.56 


0.05 


9.58 


128.0 


0.97 


3.74 


1.90 


2.05 


0.08 


1.54 


0.04 


2.73 


43.0 


1.60 


2.42 


0.93 


0.90 


0.03 


1.57 


0.03 


0.91 


11.6 


2.99 


1.08 


0.27 


0.25 


0.03 


1.31 


0.06 


0.22 


1.4 


4.00 


0.34 


0.07 


0.25 


0.03 


0.87 


0.06 


0.11 


0.7 



'^The correlation length, in physical units, computed by solving the equation ^(ro) = 1.0. The 
log^ was simply interpolated linearly between the radii r' and r", where ^(r') < 1 and ^(r") > 1. 

^ The correlation length, in physical units, computed by fitting the log^ vs log r through least 
squares. 
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Table 2. Evolutionary Parameters of ^hh 



z 


Xo 








7 


0"7 ^hh 


(1 Mpc) 




4^(0.2 Mpc) 
















,Ao)=(1.0,0.0) 








0.00 


4.43 


4.43 


4.42 


0.08 


2.41 


0.04 


35.70 


1.34 


188.8 


23.4 


0.52 


3.03 


2.00 


2.10 


0.07 


2.25 


0.06 


5.53 


0.24 


84.1 


7.8 


0.91 


2.67 


1.40 


1.41 


0.06 


2.16 


0.06 


1.91 


0.09 


47.5 


2.9 


1.48 


2.44 


0.99 


0.93 


0.04 


2.17 


0.07 


0.93 


0.07 


24.4 


1.4 


2.85 


2.24 


0.58 


0.55 


0.04 


2.26 


0.14 


0.33 


0.04 


6.6 


0.6 


4.87 


2.48 


0.42 


0.40 


0.04 


2.29 


0.14 


0.15 


0.03 


4.0 


1.0 














,Ao)=(0.2,0.0) 








0.00 


5.14 


5.16 


5.02 


0.08 


1.93 


0.02 


21.22 


0.79 


113.3 


9.7 


0.50 


4.41 


2.94 


2.99 


0.07 


1.79 


0.03 


7.44 


0.35 


61.2 


2.7 


0.99 


4.29 


2.15 


2.13 


0.06 


1.71 


0.04 


3.38 


0.15 


43.7 


3.5 


1.40 


4.14 


1.73 


1.72 


0.04 


1.66 


0.03 


2.25 


0.08 


34.8 


3.2 


2.74 


4.14 


1.11 


1.08 


0.02 


1.60 


0.02 


1.13 


0.06 


17.3 


1.4 


5.00 


4.45 


0.74 


0.73 


0.01 


1.60 


0.02 


0.61 


0.04 


7.6 


0.7 














,Ao)=(0.2,O.J 










0.00 


6.96 


6.96 


7.08 


0.10 


1.99 


0.02 


41.81 


2.59 


225.9 


26.3 


0.49 


7.54 


5.06 


4.73 


0.20 


1.86 


0.06 


19.87 


1.34 


175.4 


16.8 


0.97 


7.02 


3.56 


3.53 


0.09 


1.67 


0.03 


7.90 


0.27 


120.7 


8.5 


1.60 


6.79 


2.61 


2.54 


0.06 


1.49 


0.03 


3.74 


0.20 


66.2 


7.6 


2.99 


8.58 


2.15 


2.08 


0.05 


1.52 


0.03 


3.05 


0.26 


46.2 


11.9 


4.00 


11.87 


2.38 


2.21 


0.14 


1.55 


0.08 


3.32 


1.55 


43.4 


65.8 
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Table 3. Evolution of v^i^ and cJi, 



z 


H{z)r 










= (1.0.0.0) 




0.00 


100 


157 


1485 


0.52 


187 


384 


1192 


0.91 


263 


450 


903 








688 


2.85 


756 


172 


467 


4.87 


1421 


71 


344 




(J7o,Ao) 


=(0.2,0.0) 




0.00 


100 


154 


981 


0.50 


158 


265 


879 


0.99 


218 


288 


782 


1.40 


271 


347 


649 


2.74 


466 


191 


385 


5.00 


848 


107 


270 




(rio,Ao) 


=(0.2,0.8) 




0.00 


70 


140 


728 


0.49 


85 


140 


568 


0.97 


107 


211 


399 


1.60 


146 


121 


219 


2.99 


257 


57 


120 


4.00 


355 


39 


97 


^The 


velocities are 


measured at 


a physical 



constant separation of 1 Mpc in physical units, 
km sec~^. 
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Table 4. Evolution of < > and < ai^hh > 



z 


H{z)r 


< v\\ > cr< 


?Jj|> < (^l,hh > 


^«^i,hh> 






(fJoj-^o)— 


(1.0,0.0) 






0.00 


100 


116 


31 


861 


29 


0.52 


187 


318 


55 


924 


35 


1.48 


390 


355 


21 


557 


14 


2.85 


756 


241 


14 


413 


8 


4.87 


1421 


176 


14 


354 


11 






(J7o,Ao)= 


<0.2,0.0) 






0.00 


100 


110 


10 


460 


10 


0.52 


158 


243 


17 


504 


14 


1.48 


271 


268 


12 


425 


16 


2.85 


466 


226 


9 


335 


9 


4.87 


848 


196 


5 


278 


4 








<0.2,0.8) 






0.00 


70 


77 


9 


332 


11 


0.49 


84 


137 


14 


370 


9 


0.97 


107 


251 


13 


341 


10 


1.60 


146 


162 


10 


200 


8 


2.99 


257 


156 


9 


149 


9 


4.00 


355 


142 


20 


127 


11 



Note. — As in Table 5a the velocities here are also measured 
at 1 Mpc 
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Fig. 1. — The evolution of the two-point correlation function of the density field in physical 
coordinates measured at different epochs, from z = 0.0 to z ~ 5.0 (from top to bottom), for 
the = 1.0 scenario. The line is a power-law with exponent -1.8. 

Fig. 2. — The same as Figure 1 but for Clo = 0.2. 

Fig. 3. — The same as Figure 1 but for = 0.2. 

Fig. 4. — The evohition of from z = to z = measured by the €pp parameter is plotted as a 
function of Ze- Different types of lines denote different models: solid line for (1.0,0.0), dot-dashed 
line for (0.2,0.8), and dashed line for (0.2,0.0). Open and filled squares measure evolution at 1 Mpc 
and 0.2 Mpc, respectively. 

Fig. 5. — The evolution of the two-point correlation function of the mass density field (lines) and 
halo population (symbols) (^pp and ^hh respectively) in comoving coordinates for (f2o,Ao)=(l. 0,0.0). 
Epochs are marked as in Fig. 1. 

Fig. 6. — same as Figure 4 for (rio5Ao)=(0.2,0.0). 

Fig. 7. — same as Figure 4 and 5 for (i7o,Ao)=(0.2,0.8). In (a) ^^^^ was found by using friends of 
friends as the group finding algorithm and (b) ^^h as found by the spherical overdensity algorithm. 

Fig. 8. — Evolution of the two-point correlation function at 1 Mpc and 0.2 Mpc of the halo 
population (squares, Z = 0.1) and the mass density field (asterisks). 

Fig. 9. — Evolution of the halo two-point correlation function as measured by Shh- Symbols and 
lines are the same as in Figure 4. 
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